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Vortex structures in mesoscopic cylinder placed in exter- 
nal magnetic field are studied under the general de Gennes 
boundary condition for the order parameter corresponding to 
the suppression of surface superconductivity. The Ginzburg- 
Landau equations are solved based on trial functions for the 
order parameter for vortex-free, single-vortex, multivortex, 
and giant vortex phases. The equilibrium vortex diagrams in 
the plane of external field and cylinder radius and magneti- 
zation curves are calculated at different values of de Gennes 
"extrapolation length" characterizing the boundary condition 
for the order parameter. The comparison of the obtained vari- 
ational results with some available exact solutions shows good 
accuracy of our approach. 

PACS numbers: 74.60. Ec, 74.24. Ha 



I. INTRODUCTION 

Recent achievements in electronic device miniaturiza- 
tion allow one to study the mesoscopic superconducting 
samples with sizes of the order of the coherence length 
^(T). Such structures attract a considerable current in- 
terest as a possible basis for low temperature electronics. 
The superconducting state was studied experimentally 
for different-shaped samples: discs, loops, double loops, 
dots etc. It was shown that the sample shape and 

sizes affect significantly the phase diagrams of the meso- 
scopic superconductors. 

The vortex phases in mesoscopic superconductors are 
commonly studied within the framework of the Ginzburg- 
Landau theory As it is well-known from the mi- 
croscopic theory, the Ginzburg-Landau approach gives 
accurate results provided that the order parameter un- 
dergoes only slight spatial variation on the lenght scale 
of C(0). This means that the Ginzburg-Landau theory 
can be used in the temperature range not far from Tc- 
However, it is known from the experience that it is also 
able to give reasonable results beyond this limit. The 
Ginzburg-Landau solutions for axially symmetric meso- 
scopic samples (cylinders, discs) can be subdivided into 
two different types |p|-pT[|. In the first case the modulus 
of the local order parameter is axially symmetric inside 
the sample. The superconducting vortex-free state, the 
single-vortex state, and the giant-vortex state belong to 
this type of solutions. In the second case the axial sym- 
metry is broken and a vortex cluster is formed inside the 



sample (multivortex phase). This state usually appears 
at lower fields and larger sample sizes as compared to 
the giant- vortex phase 0,||,|o|jll| . Note that multivortex 
state corresponds to the Abrikosov fiux-line lattice for 
the bulk superconductors. 

The phase diagram of mesoscopic superconductor is 
strongly infiuenced by the boundary condition for the 
order parameter. In general case it is given by the de 
Gennes boundary condition p^Jl9[ |: 



n(-iV - A)tjj 



(1) 



where n is the unit vector normal to the sample sur- 
face, b is the de Gennes "extrapolation length", A is 
the vector potential, and ip is the order parameter. Here 
and below the following dimensionless variables are used: 
distances, magnetic field, and the order parameter are 
measured in units of coherence length ^(T), bulk upper 
critical field Hc2, and y^—a/P, respectively, with a, (3 
being the Ginzburg-Landau coefficients. The "extrapo- 
lation length" b has a physical meaning of a length scale 
of the order parameter variation at the sample surface. 
Microscopic considerations show that b depends on the 
properties of interface, it is maximum for an ideal sur- 
face with the mirror reflection of quasi-particles and min- 
imum for the rough surface with the diffusive reflection 
|l8| For the superconductor-dielectric (or vacuum) 
interface we have & — > 00 in the former case and b ^ ^(0) 
in the latter case. The value of b depends also on the sur- 
face orientation provided that a superconducting gap is 
anisotropic. It follows from Eq. (1) that the order param- 
eter is suppressed in the vicinity of the sample surface at 
6 ~ C(0)- For the superconductor-normal metal interface 
b is always small, b ^ £,{0), because of diffusion of nor- 
mal electrons from the metal to the superconductor. The 
" extrapolation length" b in this case is a function of tem- 
perature and various characteristics of metal and inter- 
face [p"9|-^3| . There are also possibilities for the enhance- 
ment of the order parameter at the interface that can 
be described by negative b values. It can be realized by 
choosing the suitable material as a surrounding medium 
p^ , p5t , i.e., a superconductor having a higher transition 
temperature then the material of the mesoscopic sam- 
ple. Another possibility is to use a semiconductor as 
a surrounding medium, such that there is a overlap of 
the band gap of the semiconductor with the supercon- 
ducting gap. For the case of isotropic superconductor- 
vacuum ideal interface (6 — > 00), the magnetic prop- 
erties of mesoscopic cylinders and discs, their equilib- 
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rium and non-equilibrium phase diagrams were studied 
in numerous papers, see e.g. Refs. ISHlq], using differ- 
ent approaches. The vortex structures allowing for the 
enhanced surface superconductivity were studied in Refs. 
p^ , p5t . In Ref. |24| the case of mesoscopic discs was con- 
sidered within the lowest Landau level approximation, 
which first was proposed in Refs. In Ref. the 

Ginzburg-Landau equations were solved numerically and 
self-consistently for superconducting state in long cylin- 
ders. 

In this paper we focus on magnetic properties of meso- 
scopic cylinders under the general boundary condition for 
the order parameter corresponding to the opposite case 
of suppressed surface superconductivity. For this pur- 
pose, we propose a variational approach and solve the 
Ginzburg-Landau equations without straightforward in- 
tegration using trial functions for the order parameter. 
These trial functions involve a set of variational param- 
eters yielding, as we show, more accurate quantitative 
description of the spatial distribution of the order pa- 
rameter than frequently used lowest Landau level approx- 
imation. The approach is applicable to all vortex phases 
(the vortex-free, the single-vortex, the multivortex, and 
the giant vortex states) and any values of de Gennes " ex- 
trapolation length" . The comparison of variational cal- 
culations with some available exact results demonstrates 
good accuracy of our approximation. The model enables 
us to calculate the equilibrium H^-R diagram of the cylin- 
der, where i/g is the external field and R is the cylinder 
radius. The magnetization curves of the cylinder are cal- 
culated. 



II. MODEL 

Let us consider a cylindrical type-II superconductor 
placed in the uniform external magnetic field parallel 
to the cylinder axis. The sample is assumed to be much 
longer than London penetration depth A(r). Therefore, 
both the order parameter and the magnetic field are con- 
stant along cylinder axis. We use the cylindrical coordi- 
nate system with coordinates r, z and unit vectors e^. 

The system of Ginzburg-Landau equations is given by 

pel: 



H{R) = H,. 
Next, we expand all variables in powers of k: 



(4) 



n=0 



n=0 



(5) 



We substitute expansions (5) to Ginzburg-Landau 
equations (2), (3) and to boundary conditions (1), (4) 
and equate powers of k. It is easy to show that the vec- 
tor potential and the magnetic field at leading order are 
given by: 



Ao = e, 



He 



Hn — H„ 



(6) 



The order parameter at leading order ip^ is determined by 
the first Ginzburg-Landau equation (2) and the boundary 
condition (1) at H = Hg. This condition is now given 
by: 



_|_ ^ Q 

dr b 

In the next order one has from Eq. (3): 



Ho = - 



dre^ 



|V'o|'Ao + -(?/'oVV'o 



(7) 



(8) 



Thus, at leading order the magnetic field is uniform in- 
side the sample and the magnetization equals zero. Phys- 
ically, this implies that the additional field generated by 
the Meissner current and by the vortices is of the order 
of I/k^ in comparison with the external uniform field. 
Below, we shall calculate the energy of the sample at 
leading order, whereas the magnetization will be found 
at the next order using Eq. (8). It was shown in Ref. p^ ] 
that this approximation is accurate not only for high-K 
materials but also for moderate-K superconductors (e.g. 
K « 4) with sizes comparable to ^ (T) . 
We present Vo as a Fourier series: 



tpo{r,'p)^ fL,{r)exp{-iLj>f). 



(9) 



|V'|^V'-V' + (»v 



rot H =— AlV'P 



(2) 



(3) 



where H and ip are the dimensionless local magnetic field 
and the order parameter (H = rot A, H = He^); k = 
A(T)/^(T) is the Ginzburg-Landau parameter. Equa- 
tions (2) and (3) must be supplemented by the boundary 
conditions for the order parameter (1) and the magnetic 
field: 



For the axial symmetric distribution of the modulus of 
the order parameter inside the sample the only one term 
in Eq. (9) is nonzero. The vortex-free, single-vortex, and 
giant vortex phases with angular quantum momentum L 
correspond to the harmonics with Lj = 0, 1, and L, re- 
spectively. The modulus of the order parameter in the 
multivortex phase is not axially symmetric. Note that 
in this case the symmetry of the vortex configuration 
imposes a restriction on functions /l, ■ some of these 
functions equal zero. It was shown in Refs. that 
taking into account only two main terms in right-hand 
side of Eq. (9) is enough for an accurate calculation of 
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the energy of the mesoscopic superconductor in multivor- 
tex state. The vortex cluster with L vortices on a ring 
and no vortex at the axis (ring-like configuration) can be 
described as a mixture of two components with Li — Q 
and L2 = L ((0 : i) state). The vortex cluster with one 
vortex at the cylinder axis and (L — 1) vortices on a ring 
corresponds to Li = 1 and L2 = L ((1 : L) state). The 
contribution from other harmonics is small and can be 
neglected, when we consider few-fluxoid cylinders that 
can accommodate just few vortices 

Using Eqs. (2) and (9) it can be easily shown that each 
function f^. (r) has the following asymptotic at r — + 0: 



(10) 



Besides, each function f^. (r) meets boundary condition 
(7). These conditions for Jl [r] are valid both for the 
giant vortex and the multivortex phases. It is a rather 
complicated task to find /l (r) from the first Ginzburg- 
Landau equation (2) due to its non-linearity. Instead of 
the straightforward integration of this equation, it is pos- 
sible to use trial functions for the coordinate dependence 
/l (r)- Notice that different variational procedures al- 
lowing one to solve approximately the Ginzburg-Landau 
equations were used in numerous papers for mesoscopic 
@,|0|jl|, bulk |2|-|3ltl, and different-shaped || super- 
conductors. One can easily show that if Eq. (9) is a solu- 
tion of the first Ginzburg-Landau equation, each function 
/l . (r) can be represented as the following series expan- 
sion in powers oir/R without loss of generality: 



exp ~qL^ 



i?2 



00 
ti=0 



r y 

r) 



V2u 



(11) 



where are constants. The value of qL can be found 
from boundary condition (7): 



9i, = 771 + o 



R_ 

2b 



2 



2u) 



y 



qPu 



(12) 



Our approach is to consider the coefficients pu' as vari- 
ational parameters and to minimize the free energy with 
respect to pt^ . The exponential prefactor in Eq. (11) 
takes into account the suppression of the order param- 
eter at the contact with a surrounding material. First 
term in the expansion describes the behavior of the or- 
der parameter in the central part of the sample. Other 
terms specify the order parameter in the sample as a 
whole. As it will be shown below, taking into account 
only first three terms in Eq. (11) is enough for the ac- 
curate calculation of the order parameter distribution, 
and we will use this approximation in all calculations of 
magnetization and H^-R diagrams. Note that trial func- 
tion (11) was used in Ref. |32| for the study of surface 
superconductivity in samples of different complex shapes 
placed in vacuum (6 — > 00). 

The Ginzburg-Landau functional for the Gibbs free en- 
ergy G of the cylinder can be written as a sum of two 



contributions, Gb and Gs- The former is the bulk energy 
of the sample and the latter is the surface energy. These 
contributions are given by pq,E3^,E6| : 



Gh 



1 



-iV- A)V|^ 



(13) 



(14) 



The integration in Eqs. (13) and (14) is performed over 
the sample bulk and surface, respectively. Note that the 
general boundary condition for the order parameter (1) 
can be obtained phenomenologically by minimization of 
the free energy functional G — Gt+Gg with respect to the 
order parameter tp and the vector potential A ||2^,^,^ . 

Substituting expansion (9) to Eqs. (13) and (14) and 
taking into account Eq. (6) we obtain the expression for 
the energy of the multivortex state (per unit length of 
the cylinder): 



Gh = 27r 



rdr 



/El /li 



dr 



dfL2 
dr 



{R)). 



(15) 



(16) 



The energy of any state having axially symmetric mod- 
ulus of the order parameter with angular momentum L 
reduces to Eqs. (15) and (16). In this case we must put 
Li = L, = 0. Using Eqs. (11) and (12) we find 
the energy G from Eqs. (15) and (16) by a straight- 
forward integration as an explicit function of variational 
parameters P2^ and Pq^ , p\^ P2^ ■ The result- 

ing expression, however, is rather cumbersome and we 
do not present it here. Finally, values of the variational 
parameters at each R are found numerically by the min- 
imization of the free energy. This procedure yields the 
local order parameter and the energy of the cylinder. 

Knowing the local order parameter we can calculate 
the magnetization. It is given by: 



AttM ^<H > -He 



(17) 



where < H > is the averaged magnetic field over the 
superconductor volume. Taking into account Eq. (8) 
and expansion (9) we obtain: 
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At:M = 



rdr / dx 



fUx) 



HpX L', 



(18) 



Functions /lj and /^^ are found by the method described 
above. Hence, one can calculate the magnetization using 
Eq. (18). In the following section we apply the developed 
approach for the analysis of the behavior of the cylinder 
in the external field. 



III. RESULTS AND DISCUSSION 

Comparing the energies of different states one can cal- 
culate the equilibrium H^-R diagram of the cylinder. The 
results of our calculations are shown in Fig. 1 for differ- 
ent h values: h = 1 (a), & = 2.5 (b), 6 = 5 (c), and 
& ^ oo (d). The latter case corresponds to the isotropic 
superconductor-vacuum ideal interface, and was studied 
in Refs. p|,p|,p^,pT jl3 - 15[| . Curve 1 shows the transition 
from the normal to the superconducting state (the sur- 
face critical field -ffcs)- The oscillatory behavior of the 
function Hcz{R) is caused by the fact that the transi- 
tion occurs from the normal to the giant vortex states 
with different angular quantum moments L depending 
on the cylinder radius. Besides, the function Hc3{R) de- 
pends appreciably on the value of b: the value of HcsiR) 
decreases with decrease of 6. At i? ^ oo the dependence 
HcsiR) tends to the surface critical field for the half-space 
sample, which was calculated in Ref. as a function 
of b. Note that Hc3 (oo) = 1.695 for 6 — > oo and Hcsioo) 
= 1 for 6 ^ |l|,§. 

Below Hc3 the transitions between different giant vor- 
tex states take place. Solid lines show the phase bound- 
aries between the states with different vorticities that 
are the sums of angular quantum moments of all vortices 
and giant vortices. In the giant vortex state the order 
parameter is strongly suppressed in the inner part of the 
cylinder, and this state can be referred to a surface su- 
perconductivity. For illustration, the spatial distribution 
of the order parameter in the giant vortex phase with 
L = 2 is plotted in Fig. 2 at 5 = 1, i? = 3.9, He = 0.9 
(a) and b oo, R = 4.56, — 0.5 (b) (solid lines). In 
the former case the order parameter is also suppressed 
at the sample surface because of the small b value (e.g., 
superconductor-normal metal interface). To check the 
accuracy of our approach we took into account next sev- 
eral terms in expansion (11), thus increasing the number 
of variational parameters. As we found, this practically 
did not change the calculated order parameter for almost 
all points of Hf.-R diagram shown in Fig. 1. This result 
implies that our variational calculations are close to the 
exact solutions of the Ginzburg-Landau equations since 
expansion (11) is written without any loss of generality. 
Therefore, one can find Ginzburg-Landau solution with 



any desired accuracy (for the states with axially symmet- 
ric modulus of the order parameter) allowing for enough 

number of variational parameters pu' . For multivortex 
states the accuracy is limited by the fact that we take into 
account only two main harmonics in the Fourier expan- 
sion (9). In Fig. 2 we also plotted the spatial variation of 
the order parameter calculated within the lowest Landau 
level approximation (dashed lines). In this approxima- 
tion the order parameter is proportional to the eigen- 
function of the kinetic energy operator corresponding to 
the lowest eigenvalue. We found that this approach re- 
mains very accurate not far from Hc2 and i/cs (see Fig. 
2(a)). At lower fields the results of the lowest Landau 
level approximation are not so accurate (Fig. 2(b)). 

As follows from Fig. 1, the superconducting state does 
not nucleate at very small cylinder radii, smaller than 
some critical radius, and the sample is in normal state 
at any applied field. The critical radius tends to zero at 
6 — > oo. There is also the interval of R for each 5, when 
the vortex phase does not nucleate, and the transition oc- 
curs from the normal to the superconducting vortex-free 
state. Every vortex phase with vorticity L > 1 can exist 
in the form of the giant vortex or the multivortex config- 
uration. The dashed lines on Fig. 1 show the boundaries 
between these states. Below these lines for given L > 1 
the multivortex state has the lowest energy, and above 
these curves the giant vortex state becomes more ener- 
getically favorable. In equilibrium state the multivortex 
phase can exist if the applied field is smaller than 1 (i/c2 
in dimensional units) and if the radius of the cylinder is 
large enough. For each b there exists an interval of small 
cylinder radii, when the multivortex phase is energeti- 
cally unfavorable as compared to the states with axial 
symmetric distribution of the modulus of the order pa- 
rameter. Curve 2 in Fig. 1 shows the lower critical field 
that corresponds to the equilibrium boundary between 
the vortex-free and the single-vortex states. 

With increasing the external field the cylinder can fol- 
low rather complex set of phase transitions. It can come 
from the giant vortex to the multivortex states and then 
back to the giant vortex phase. For example, at i? = 4.4, 
6=1 the following set of transitions occurs: 0^1^ 
(0:2) 2 ^ (0:3) ^ 3 (0:4) -^4. At i? = 4.1, 6 = oo 
the set of transitions is: — > 1 — > 2 — > (0:3) — > 3 
(0:4) ->4->5-^6-^7^8^9^10^11. 

The transitions between the phases with different vor- 
ticities are always discontinuous, they occur when the 
energies of different states become equal. The transi- 
tions between the states with the same vorticity may be 
continuous as well as discontinuous. Continuous phase 
transitions occur between the multivortex and the gi- 
ant vortex states. In this case, with increasing of the 
applied field the intervortex distances decrease, and vor- 
tices merge into the giant vortex located at the cylinder 
axis. The similar results were obtained in Refs. [7|-[Tl|] 
for the case 5 = oo and in Ref. at 6 < (enhanced 
surface superconductivity). In our calculations we also 
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took into account a possibility of existence of the mul- 
tivortex clusters with central vortex. It turned out that 
such configurations can be more favorable energetically 
than giant vortex states in some regions of H^-R diagram 
shown in Fig. 1. However, in all these cases the ring- like 
clusters have the lower energy. For cylinders thicker than 
shown in Fig. 1 we found that at any 5 > the ground 
state can be represented by configurations with the cen- 
tral vortex. The transitions between multivortex states 
of the same vorticity with and without central vortex are 
discontinuous. 

With increasing of the cylinder radius the number of 
vortices, which it is able to accomodate, increases, and 
finally vortex array transforms to the classical triangular 
Abrikosov flux-line lattice far from the surface. The first 
step on this way is the appearance of the clusters with 
central vortex. However, we do not analyze here the tran- 
sition from the mesoscopic to the macroscopic behavior 
and restrict ourselves on few-fluxoid cylinders, which can 
accomodate only few vortices before the thansition to the 
normal state and whose magnetic properties can be de- 
scribed in terms of mixture of only two harmonics in Eq. 
(9). 

Now we find the magnetization of the cylinder using 
Eq. (18). The results of our calculations at k = 5 are 
presented in Fig. 3(a) for b — 1, R = 4.625 and in Fig. 
4(a) for b oo, R = 4.05. In these cases the cylinders 
accommodates the giant vortices with maximum angular 
quantum moments equal to 5 and 11, respectively, before 
the transition to the normal state. Jumps in the magne- 
tization correspond to the transitions between the states 
with different vorticity. It is interesting that in the state 
with suppressed surface superconductivity (Fig. 3(a)) 
the discontinuity of the field dependence of the magne- 
tization is less pronounced, especially at high fields. In 
the first case (Fig. 4(a)) at Hf. « 0.61 (0.61i?c2 in dimen- 
sional units) the transition occurs from the multivortex 
state with 2 vortices to the giant vortex state with angu- 
lar quantum momentum 2 (2 — > (0:2)). In the second case 
(Fig. 4(b)) the transition 3 (0:3) occurs at « 0.74 
and the transition 4 (0:4) occurs at « 0.87. All 
these transitions are followed by weak jumps in the slope 
of the magnetization. The behavior of the magnetization 
near the transitions is shown in Fig. 3(b) and 4(b). Solid 
lines denote the equilibrium magnetization, dashed lines 
denote the metastable magnetization corresponding to 
the giant vortex phase. The similar behavior of the mag- 
netization in the vicinity of the multivortex-giant vortex 
transitions was reported in Ref. [34] for the case b ^ oo. 

Next we discuss the accuracy of our variational pro- 
cedure. First, we compare the variational results with 
known exact solutions for the surface critical field. The 
phases with axial symmetric distribution of the modulus 
of the order parameter are always energetically more fa- 
vorable with respect to the multivortex state at applied 
field higher than the bulk upper critical field (sec phase 
diagrams on Fig. 1). In the vicinity of the surface critical 
field the first Ginzburg-Landau equation (1) can be lin- 



earized. A resulting equation has the following analytical 
solution (s): 

/.(i?) exp j $ L + 1, ^) , 

(19) 

where $ is Kummer function. Function (19) must meet 
the boundary condition (7). This yields the transcenden- 
tal equation for the surface critical field allowing one to 
find HcsiR) exactly. The comparison of the variationally 
calculated Hcs with this exact dependence shows good 
agreement with an accuracy better than 1% for all values 
of b and R under study. The lower critical field of the 
cylinder Hd versus R was calculated numerically in Ref. 
]l2[ at b —> 00. The comparison of this result with our 
dependence Hci{R) reveals the same accuracy. Thus, our 
results appear to be a good approximation to the exact 
Ginzburg-Landau solutions for the mesoscopic cylinders. 

In summary, we analyzed the superconducting state in 
long mesoscopic cylinder with suppressed surface super- 
conductivity. An asymptotic expansion was used to sim- 
plify the Ginzburg-Landau equations at high and moder- 
ate values of k, and the simplified equations were solved 
by variational method. The equilibrium H^-R diagram 
of the cylinder were obtained, where is the external 
field and R is the cylinder radius, at different values of 
" extrapolation length" . We showed that magnetic prop- 
erties of the cylinder depend appreciably on the value of 
" extrapolation length" and studied the evolution of He-R 
diagram with changing of " extrapolation length" . 
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FIGURE CAPTIONS 



Fig.l. Equilibrium H^-R diagram of the cylinder in the external magnetic field at 6 = 1 (a), h = 2.5 (b), 6 = 5 (c), 
6 — > cx) (d). Solid lines show the boundaries between the states with difi'erent vorticity. Dashed lines correspond to 
the boundaries between the multivortex and the giant vortex phases. Curves 1 and 2 show the surface and the first 
critical fields, respectively. The dot line denotes the bulk upper critical field. 

Fig 2. The spatial distribution of the dimensionless modulus of the order parameter inside the mesoscopic cylinder 

in giant vortex state with angular quantum momentum L = 2iith=l, R = 3.9, He = 0.9 (a) and b oo, R = 4.56, 
He — 0.5 (b) Solid lines correspond to our variational result, dashed lines to the results of the lowest Landau level 
approximation. The distance from the cylinder axis r is measured in units of the coherence length ^(T). 

Fig. 3. The equilibrium magnetization of the cylinder with radius R = 4.625 versus applied field at 6 = 1, k = 5. 
Jumps in the magnetization in Fig. 3(a) correspond to the transitions between the states with different vorticity. Fig. 
3(b) shows the behavior of the magnetization in the vicinity of the continuous phase transition at Hg « 0.61iJc2 from 
the multivortex state with 2 vortices to the giant vortex phase with vorticity L = 2. Solid lines correspond to the 
equilibrium magnetization, dashed line shows the magnetization of the metastable giant vortex phase. 

Fig. 4. The equilibrium field dependence of the magnetization of the cylinder with radius 4.05 at 6 — » oo, k = 5. 
In Fig. 4(b) the magnetization is plotted versus applied field near the transitions from the multivortex states to the 
giant vortex phases. Solid lines correspond to the equilibrium magnetization, dashed line shows the magnetization of 
the metastable giant vortex phase. 
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Fig. 4(a) 
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